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I.  INTRODUCTION 


This  final  report  contains  the  results  of  the  project  SPC-004064  and  concentrates  mainly 
latest  results  of  measurement  of  electrical  conductivity  in  liquid  CdTe  and  under  a  defined  Cd 
overpressure.  Attention  is  also  paid  to  evaluation  of  high  temperature  galvanomagnetic  mea¬ 
surements  in  order  to  find  optimal  cooling  conditions  for  fabrication  of  near-stoichiometric 
single  crystals  from  the  point  of  view  of  chemical  self-diffusion  of  Cd.  Due  to  the  fact,  that 
interest  in  fabrication  of  large,  homogeneous  high-resistivity  CdTe  and  (CdZn)Te  for  gamma 
ray  detectors  applications  strongly  increases,  conditions  for  preparation  of  such  a  material 
with  a  minimized  total  concentrations  of  defects  and  thus  with  improved  detection  prop¬ 
erties  are  looked  for.  A  possibility  to  use  measurements  of  optical  absorption  to  evaluate 
concentrations  of  shallow  absorption  levels  are  discussed  as  well.  The  topics,  which  were  in 
detail  described  in  the  Progress  report  to  this  project  are  not  repeated  in  the  Final  report. 

II.  HIGH  TEMPERATURE  ELECTRICAL  PROPERTIES  OF  LIQUID  CDTE 

A.  Introduction 

CdTe  and  (CdZn)Te  with  a  small  Zn  content  have  considerable  application  in  optoelec¬ 
tronics  devices.  These  materials  are  used  as  substrates  for  epitaxial  growth  of  HgCdTe  layers 
for  infrared  detectors,  for  a  fabrication  of  gamma  ray  detectors  and  electrooptical  modu¬ 
lators.  All  these  applications  require  large  bulk  single  crystals  with  defined  electrical  and 
optical  properties. 

The  growth  of  single  crystals  from  the  melt  represents  the  dominating  method  of  prepa¬ 
ration  of  crystals  to  the  above  mentioned  purposes.  The  most  commonly  used  experimental 
arrangement  is  the  vertical  directional  cooling  of  the  stoichiometric  melt.  The  routine  fab¬ 
rication  of  CdTe  and  (CdZn)Te  crystals  with  pre-defined  electrical  and  optical  parameters 
has  not  been  fully  managed.  This  fact  is  caused  mainly  by  relatively  non-standard  physic¬ 
ochemical  properties,  such  as  a  high  partial  pressure  of  a  more  volatile  component  around 
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the  melting  point,  conservation  of  bonds  between  Cd  and  Te  atoms  in  liquid  CdTe  and  a 
low  thermal  conductivity,  which  substantially  differs  in  the  liquid  and  in  the  solid  state. 
Especially  the  first  two  properties  influence  the  growth  of  single  crystals  from  the  melt  and 
therefore  we  will  discuss  them  in  the  next  chapter  more  deeply.  The  bad  thermal  conduc¬ 
tivity  of  the  solid  results  only  in  a  gradual  slowdown  of  the  crystal  growth,  especially  in  the 
last  stage,  when  the  heat  outlet  of  the  latent  heat  of  crystallization  is  weak. 

Therefore  a  method  is  looked  for,  how  to  receive  maximum  information  about  the  whole 
process  of  preparation  of  single  crystals  from  the  melt  and  especially  about  the  way,  how 
the  growth  and  quality  of  single  crystals  is  influenced  by  external  growth  conditions  in  order 
to  optimize  the  process.  So  far  very  little  is  known  about  the  state  of  the  melt  before  the 
beginning  of  crystallization  and  about  the  respective  process  of  solidification  of  the  melt. 
Therefore  we  concentrated  our  effort  on  these  two  topics. 

It  came  out,  that  the  measurement  of  electrical  conductivity  is  suitable  for  these  pur¬ 
poses.  Originally  we  wanted  to  study  also  the  Hall  effect,  which  in  general  comprises  more 
information  about  the  current  carriers,  but  due  to  the  fact,  that  this  effect  is  orders  of  mag¬ 
nitude  smaller  than  the  conductivity,  it  is  strongly  influenced  by  the  noise,  which  dominates 
the  signal  at  temperatures  above  lOOO^C. 

B.  The  properties  of  liquid  CdTe  and  preparation  of  crystals  from  the  melt 

The  process  of  crystal  growth  from  the  melt  can  be  divided  in  three  stages. 

1st  stage  -  the  melt 

Two  phases  coexist  -  the  gas  phase  (G)  and  the  liquid  phase  (L),  when  the  state  of  the 
melt  is  deflned  by  temperature  T  and  a  pressure  of  one  of  the  components  e.g.  Cd  {Pcd)- 
The  congruent  melting  point  of  CdTe  is  Tm=1092''C.  Relatively  little  is  known  about  the 
properties  of  liquid  CdTe.  To  describe  the  liquid  a  Jordan’s  model  of  regular  associated 
solutions  is  applied.  Model  calculations  yield,  that  the  stoichiometric  liquid  is  only  partially 
dissociated  (about  5  percent)  into  Cd  and  Te  atoms  [1].  The  conservation  of  bonds  in 
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liquid  CdTe  was  confirmed  by  neutron  diffraction  experiments  [2],  Hall  effect  [3],  electrical 
conductivity  [3-6].  Also  theoretical  calculations  [7]  suppose  existence  of  Te  chains,  to  which 
Cd  is  bonded.  Muhlberg  [8]  found  out  during  the  growth,  that  a  better  quality  single  crystals 
can  be  obtained  from  highly  superheated  melts,  when  the  dissociation  of  the  liquid  is  higher. 
Most  experimental  measurements  were  so  far  performed  only  in  dependence  of  temperature, 
without  a  concrete  knowledge  of  a  partial  pressure  of  one  of  the  components.  That  is  why  the 
experimental  dependences  obtained  such  a  way  must  be  viewed  with  a  great  care.  As  will  be 
shown  later  in  the  report,  the  dependence  on  the  partial  pressure  of  one  of  the  components 
is  strong. 

2nd  stage  -  the  solidification 

This  stage  is  characterized  by  coexistence  of  three  phases  -  the  gas(G),  the  liquid  (L) 
and  the  solid  (S).  The  state  of  both  the  liquid  and  solid  phases  is  defined  by  the  partial 
pressure  of  one  of  the  components  on  the  three  phase  line.  We  found  out  [9],  that  the  Cd 
pressure  corresponding  to  the  stoichiometric  solid  is  Pc'd=l-9atm  and  to  the  stoichiometric 
liquid  Pcd= l-2atm.  It  means,  that  stoichiometric  solid  is  fabricated  by  solidification  of 
the  Cd  rich  liquid,  while  solidification  of  stoichiometric  liquid  results  in  Te  rich  solid.  The 
dominating  component  in  the  vapor  phase  is  in  both  cases  Cd. 

Another  interesting  point  on  the  three  phase  line  is  the  point,  at  which  the  gas  phase 
is  stoichiometric,  it  means  when  congruent  sublimation  occurs.  This  pressure  is  approx 
P^^=0,12atm.  At  this  pressure  at  the  given  temperature  is  the  total  pressure  above  the 
completely  solidified  liquid  minimal  and  the  material  sublimes  strongly.  Our  experiments 
have  shown,  that  e.g.  electrical  conductivity  can  be  measured  only  at  pressures  Pcd  >1*1 
atm,  when  sublimation  does  not  influence  the  performed  measurements. 

3rd  stage  -  the  solid 

Two  phase  coexist  during  the  cooling  of  the  solid  -  the  gas  (G)  and  the  solid  (S).  The 
state  of  the  solid  is  in  this  case  determined  by  temperature  and  a  pressure  of  one  of  the 
components,  e.g.  Cd.  We  studied  the  processes  occurring  within  the  range  of  coexistence  of 
S  and  L  within  the  last  years.  We  found  out  to  evaluate  from  annealing  experiments  and  from 
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in  situ  [9-11]  galvanomagnetic  measurements,  when  we  measured  electrical  conductivity  and 
Hall  effect  the  following  conclusions 

a)  to  determine  the  line  Pcd{S)  in  the  P-T  diagram,  which  corresponds  to  the  stoichio¬ 
metric  solid,  it  means  the  mode  suitable  for  cooling  of  the  crystal  to  keep  its  stoichiometry 

b)  to  clarify  the  processes  of  diffusion  which  define  the  process  of  establishing  of  the 
equilibrium  in  the  crystal,  it  means  an  optimized  time  procedure  of  crystal  cooling. 

It  is  clear,  that  information  about  the  equilibrium  of  the  liquid  and  gas  phases  are 
insufficient.  No  authors  studied  both  the  temperature  and  the  pressure  dependence  of  any 
parameter,  which  could  provide  information  about  cooling  of  the  liquid  including  the  phase 
transition  between  the  solid  and  the  liquid.  We  tried  to  study  the  electrical  conductivity  in 
dependence  of  T  and  Pcd  In  the  melt  and  during  the  solidification  of  the  melt. 

C,  Experimental  setup 

One  of  the  main  goals  of  the  present  project  is  to  measure  galvanomagnetic  properties 
in  liquid  CdTe  and  CdZnTe  under  a  defined  Cd  overpressure.  To  complete  this  difficult 
task  we  constructed  a  new  experimental  setup  based  on  establishing  a  well  defined  partial 
pressure  of  the  volatile  component  in  a  semi-open  system  originally  described  by  Van  Doom 
[12],  which  he  used  in  the  Philips  laboratories.  The  principle  of  the  method  is  apparent  from 
Fig.  1,  the  total  arrangement  of  the  experiment  from  Fig.  2. 

The  CdTe  sample  is  located  in  the  plateau  of  constant  temperature.  The  pressure  of  Cd 
vapor  is  kept  by  evaporation  of  Cd  from  the  source.  Cd  condenses  on  the  wall  at  the  place, 
where  the  temperature  reaches  the  boiling  point  of  Cd.  The  liquid  Cd  then  flows  along  the 
wall  to  places  with  a  higher  temperature  where  it  again  evaporates.  This  cycle  is  similar 
to  the  way,  how  the  diffusion  pump  works.  The  place  of  Cd  condensation  depends  on  the 
pressure  of  Ar,  which  is  in  equilibrium  with  the  Cd  pressure.  The  pressure  of  Ar  is  kept  on 
a  constant  value  with  a  pressure  gauge,  which  enables  an  automatic  control  of  pressure  with 
a  high  precision  up  to  1.6  atm.  The  detailed  arrangement  of  the  measuring  cell  is  apparent 
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from  Fig.  3. 


FIG.  1.  The  scheme  of  the  measuring  cell  within  the  Van  Doom  adjustment. 


FIG.  2.  Schematic  drawing  of  the  furnace  control.  Both  temperature  and  pressure  arc 
controlled  by  Eurotherm  2704  dual  loop,  high  accuracy  and  stability  process 
controller. 
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FIG.  3,  The  measuring  ampoule  for  the  measurements  of  conductivity  of  liquid  and  solid  CdTe 
and  (CdZn)Te  in  the  Van  der  Pauw  configuration.  Four  graphite  contacts  are  seen  beneath,  the 
sample  is  at  tho  bottom  of  the  ampoule.  The  Cd  source  is  above  the  bottleneck. 

The  measurements  were  performed  in  both  the  DC  and  the  AC  modes,  when  better 
results  were  obtained  in  the  DC  mode.  We  concentrated  on  measurements  of  the  electrical 
conductivity,  because  the  measurements  of  the  Hall  coefficient  were  affected  by  a  great 
experimental  error.  Furthermore  the  measured  dependences  showed  a  similar  slopes  and 
trends  as  the  data  from  the  measurements  of  electrical  conductivity.  It  comes  out,  that 


should  it  be  necessary  to  study  the  Hall  effect  in  the  liquid  CdTe,  a  new  type  of  the  measuring 
cell,  which  enables  to  measure  very  thin  layers  will  have  to  be  constructed.  Preliminarily 
we  start  to  work  on  this  topic. 

D.  Electrical  conductivity  of  liquid  CdTe 

The  dependence  of  electrical  conductivity  in  dependence  of  temperature  and  Cd  partial 
pressure  is  shown  on  Fig.  4.  A  fully  new  result  is  the  dependence  of  the  degree  of  supercooling 
of  the  melt  on  Cd  partial  pressure.  This  result,  which  is  apparent  from  the  figure,  is  to  some 
extent  unexpected.  The  supercooling  AT  decreases,  when  Pcd  increases.  For 
1.4,  and  1.6  atm,  AT=8,  13,  23®C,  respectively.  The  cooling/heating  rate  15^C/hour  was 
used  in  the  measurement.  We  can  compare  our  results  with  [13],  where  the  dependence  of 
supercooling  on  Te  mole  fraction  in  Te-rich  CdTe  x  is  reported.  For  the  decreasing  x=0.55, 
0.51,  0.5  the  supercoling  decreases;  AT=100,  75,  38''C,  respectively.  We  can  conclude,  that 
both  in  Te-rich  and  in  Cd-rich  material  AT  decreases,  when  the  content  of  Cd  in  CdTe 
increases.  This  fact  has  an  impact  to  the  crystal  growth  because  the  large-diameter  single 
crystals  can  be  grown  only  in  the  regime  of  minimized  supercooling.  This  condition  is  better 
fulfilled  in  the  Cd-rich  material,  than  in  the  stoichiometric  one.  The  crystal  growth  at  high 
pressure  of  Cd  should  be  thus  examined  to  verify  such  a  possibility.  In  case  of  stoichiometric 
material  the  supercooling  must  be  minimized  decreasing  the  cooling  rate  or  by  a  proper 
arrangement  of  the  crystal  growth  apparatus. 
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FIG.  4.  The  conductivity  of  solid  and  liquid  CdTe  vs  temperature  at  various  Cd  pressures 
=1.3,  1.4,  1,5  atm.  The  direction  of  the  cooling/heating  loop  is  indicated  by  arrows,  cl  and 
as  show  fit  of  a  in  the  liquid  and  solid,  ateor  is  theoretical  conductivity  of  solid.  The  melting  point 
Tm  is  labelled  as  well. 

The  guide  full  lines  in  Fig.  4  describe  a  of  solid  and  liquid.  In  case  of  as  in  the  solid, 
which  is  defined  by 

as  —  3.3  X  10^exp(-1.2eV/A;6T)  (fi'^cm”^),  (1) 

the  interpretation  is  not  clear.  The  comparable  theoretical  curve  ateor  calculated  within  our 
model  is  higher  and  its  slope  is  less,  only  0.8  eV.  The  conductivity  in  the  liquid 

aL  =  9.1  X  I0^exp(-0.8eV/fc^,T)  (Q-^cm"^)  (2) 

agrees  well  with  [6],  where  the  slope  0.85  eV  was  reported.  The  ai  can  be  properly  ap- 
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proximated  within  our  model  of  native  defects  in  the  solid  [3],  if  we  set  deep  donor  level 
with  density  ~  10^^  cm”^  about  1  eV  below  the  bottom  of  the  conduction  band.  On  such 
a  level  Te  antisite  Tecd  is  often  assumed.  Concluding  from  that,  the  liquid  CdTe  behaves 
electrically  as  semiconductor  with  high  density  of  deep  defects,  mostly  donors.  These  defects 
are  probably  dissociated  free  Te  or  Cd  atoms.  Their  energy  levels  are  deep  below  the  Fermi 
level,  which  results  about  200  meV  below  the  conduction  band.  The  transport  properties 
are  dominated  by  electrons  in  conduction  band  excited  from  the  donors.  Increasing  temper¬ 
ature,  the  fraction  of  dissociated  species  increases  and  the  material  is  expected  to  become 
a  metal  -  second  order  phase  transition.  Looking  for  such  a  crossing  point  (semiconductor- 
semimetal)  is  desirable  to  understand  the  statistics  of  the  defects  and  transport  mechanisms 
in  the  liquid. 

A  comparison  of  our  experiment  with  other  relevant  data  is  presented  in  Fig.  5.  In  case 
of  [4]  and  [5]  the  measurements  without  defined  Pcd  were  done.  These  curves  have  similar 
features  as  our  a  measured  at  the  cooling  at  1.3  atm.  As  we  have  observed,  the  sublimation 
of  the  material  is  high  at  this  pressure  and  the  phase  transition  is  difficult  to  be  proven. 
Much  better  results  are  published  in  [6],  which  presents  as  at  Pcrf=l-42  atm  very  close  to  our 
theoretical  value  and  phase  transition  including  supercooling  well  pronounced.  The  melting 
point  is  about  10°C  below  correct  value.  Indirect  method  to  measure  electrical  conductivity 
in  closed  system  by  multifrequency  eddy  current  sensor  is  reported  in  [14].  This  method 
looks  optimum  for  the  examination  of  the  phase  transition  and  to  collect  the  essential  data 
for  the  crystal  growth.  The  absolute  value  of  a  is  affected  by  a  systematic  error  and  a  is 
less  than  data  from  direct  measurements. 
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FIG.  5.  Comparison  of  our  measurement  of  cr  at  1.4  atm  (a)  with  data  reported  by  other  workers, 
(b)  [4],  (c)  [5],  (d)  [6],  (e)  [14] 


E.  Future  plans 

The  reproducible  growth  of  high-quality  single  crystals  requires  a  detailed  understanding 
of  the  relationships  between  the  controllable  growth  parameters  and  the  properties  of  CdTe, 
(CdZn)Te,  especially  in  solid  and  liquid  near  three  phase  line  in  P-T  diagram  and  phase 
equilibria  during  the  crystal  growth.  High  temperature  in- situ  electrical  measurements  at 
defined  thermodynamic  conditions  are  an  effective  tool  to  obtain  relevant  data  from  this 
field.  Based  on  the  results  discussed  above  we  plan  the  following  experiments: 

1.  To  measure  electrical  conductivity  of  liquid  CdTe  at  high  Cd  pressures  (Pcd=h6  -  10 
atm).  The  maximum  Pea  at  the  three  phase  line  in  the  phase  diagram  will  be  looked  for 
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and  its  neighborhood  studied  in  detail.  The  properties  of  liquid  and  solid  near  this  critical 
point  will  be  evaluated. 

2.  To  study  the  influence  of  the  overheating  and  of  the  cooling  velocity  to  the  supercooling 
of  the  melt  and  to  propose  an  optimum  overheating  for  the  crystal  growth. 

3.  To  study  in  the  liquid  CdTe  the  point,  where  the  semiconducting  behaviour  changes 
to  the  semimetallic  one.  To  complete  the  P-T  phase  diagram  by  the  line  at  which  this  phase 
transition  of  second  order  occurs. 

4.  There  is  an  open  question  about  the  dependence  of  the  growth  direction,  in  our 
crystals  of  type  (111),  which  changes  to  the  direction  (100),  if  Pcd  is  increased.  We  would 
like  to  prove  and  explain  such  observation  and  to  specify  the  limiting  Pcd,  which  determines 
such  conversion. 


III.  CHEMICAL  SELF-DIFFUSION  IN  CDTE 

A.  Introduction 

The  chemical  diffusion  coefficient  D  is  a  key  parameter  characterizing  an  annealing  time, 
which  is  necessary  to  attain  the  thermodynamic  equilibrium  at  given  temperature  or  to 
produce  a  desired  gradient  in  composition.  Chemical  self-diffusion  experiments  in  near  stoi¬ 
chiometric  bulk  CdTe  are  carried  by  in-situhigh  temperature  electrical  measurements,  where 
the  response  time  of  the  free  carrier  concentration  or  electrical  conductivity  corresponding 
to  step  changes  of  ambient  partial  Cd  pressure  {Pcd)  are  measured,  for  a  review,  see  [15]. 
The  only  relevant  experimental  data  on  undoped  samples  are  found  in  [6],  which  reports  on 
the  diffusion  at  temperatures  T=700-900^C  and  in  [16]  for  T=600-800°C. 

The  aim  of  this  paper  is  to  check  the  above  mentioned  data  and  to  repeat  the  mea¬ 
surements  in  CdTe  of  a  better  purity  than  they  were  performed  before.  No  experimental 
measurements  have  been  made  since  ~  1980  and  new  experiments  are  highly  desirable  in 
this  field. 
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B.  Experiment 


Undoped  CdTe  single  crystals  were  produced  by  the  Vertical  Gradient  Freeze  method 
[17]  at  Cd  overpressure  1-1.3  atm.  6N  purity  starting  elements  were  used  in  all  growth 
runs.  In  situ  galvanomagnetic  measurements  of  conductivity  a  and  Hall  effect  by  the  Van 
der  Pauw  method  on  samples  with  dimensions  10x10x1.5  mm^  were  performed  in  two  zone 
furnace  at  temperatures  up  to  1050''C  and  Pcd  controlled  by  a  temperature  Tcd  of  the  cold 
zone  [3]. 

Time  evolution  of  samples  after  step-like  change  of  Pcd  in  the  interval  of  Pcd  one  order 
of  magnitude  below  and  above  the  ideal  stoichiometry  line  was  measured  and  the  chemical 
diffusion  coefEcient  of  cadmium  telluride  D  was  determined  as  a  function  of  temperature  in 
the  temperature  range  r=500"700°C.  The  delay  of  the  furnace  to  reach  stable  state  after 
the  step-like  change  of  Tcd  was  about  10  min. 

C.  Theory 

The  analysis  of  high  temperature  measurements  allowed  us  to  establish  the  region  in  the 
pressure-temperature  (P-T)  phase  diagram  to  obtain  intrinsic  CdTe  and  to  correlate  them 
with  the  stoichiometry  of  the  samples.  The  position  of  the  theoretical  stoichiometry  line 
(the  so  called  p-n  line)  in  the  form 

P^  =  8  X  10^exp(-1.76  x  lOV^)  (atm).  (3) 

was  found  [9].  Simultaneously,  the  galvanomagnetic  properties  of  CdTe  were  used  to  estab¬ 
lish  formation  energies  and  entropy  terms,  which  determine  the  properties  of  native  point 
defects  [3]. 

Deviating  from  native  point  defects  which  represent  excess  or  deficit  of  Cd  are  created 
and  diffuse  in  the  sample.  The  favorite  electrically  active  native  defects,  which  dominate  in 
near  stoichiometric  samples  CdTe  are  Cd  interstitial  (Cd/,  divalent  donor)  and  Cd  vacancy 
(Vcd,  divalent  acceptor)  [3,18].  In  view  of  the  fact,  that  maximum  dissolvable  excess  of 
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Cd  or  Te  in  solid  CdTe  is  about  four  times  larger  than  the  defect  density  obtained  from 
electrical  measurements  [18],  neutral  native  defects  are  also  present  in  the  sample.  Their 
concentration  depends  much  strongly  on  Pcd  than  it  is  in  case  of  charged  native  defects  and 
they  prevail  in  CdTe  in  the  regions  of  P-T  diagram  near  the  three-phase-line  [9] .  Charged 
and  neutral  defects  should  be  taken  into  account  to  describe  the  chemical  diffusion  in  CdTe 
properly.  They  can  act  both  actively  transferring  material  within  the  sample  and  passively 
as  traps  for  the  fast  moving  defects. 

Such  a  complicated  structure  of  native  defects  results  in  a  necessity  to  assume  D{A)  as 
a  function  of  the  stoichiometry  deviation  A.  The  experimental  data  show,  however,  that 
D{A)  is  practically  independent  of  A  in  samples  near  stoichiometry  [6,16].  Also  in  our 
measurements  we  were  not  able  to  prove  any  dependence  on  A  and  we  follow  the  previous 
approaches. 

In  our  arrangement  of  a  flat  sample,  the  one-dimensional  diffusion  can  be  assumed  and 
the  diffusion  profile  of  the  stoichiometry  variations  reads 

A{x,  t)  =  (Ai  -  Ao)E“  o(-l)"  erfc  )  +  erfc  ^ 

where  the  initial  and  final  stoichiometry  deviations  Aq,  Ai,  respectively,  are  calculated 
within  our  model  [3].  They  correspond  to  A  of  the  sample  in  equilibrium  at  Pcd  before 
and  after  the  step-like  change.  L  is  the  thickness  of  the  sample  {0  <  x  <  L),  homogeneous 
sample  is  assumed  at  t  =  0. 

Near  where  the  divalent  electrically  active  native  defects  dominate,  the  defect  statis¬ 
tics  according  to  [3,9]  is  used  and  profiles  of  free  carrier  densities  [electrons  n(x,t),  holes 
h{x,t)]  are  calculated  via  their  relation  to  A{x,t).  Consequently,  the  conductivity  a{t)  is 


given  as 


e 

r{t)  =  7  /  t)  +  iJ.hh{x,  t)]dx, 

Jj  Jo 


where  fJie{T),  fJh{P)  are  bulk  experimental  electron  and  hole  mobilities,  respectively.  The 
diffusion  coefficient  D  is  the  only  optimized  parameter  to  fit  a{t)  to  the  experimental  con¬ 


ductivity  data. 
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Below  600°C  the  surface  conductivity  manifests  in  the  initial  fast  relaxation  of  a.  This 
effect  must  be  separated  to  distinguish  the  bulk  diffusion  properly. 

D.  Results  and  discussion 

The  chemical  self-diffusion  was  studied  at  temperatures  SOO'^C,  GOO'^C,  and  TOO'^C.  The 
results  are  summarized  in  Table  I,  where  the  relevant  data  of  other  authors  are  included. 
The  relatively  thin  samples  (1.5  mm)  used  in  our  measurements  allowed  to  measure  the 
diflFusion  at  500°C,  where  no  data  were  reported  yet.  On  the  contrary,  fast  diffusion  at 
T  ^  700^^0,  which  resulted  in  the  response  time  comparable  with  the  characteristic  time  of 
the  furnace  (10  min),  did  not  allow  to  establish  D  definitely.  The  average  chemical  diffusion 
coefficient  deduced  from  data  at  500  and  GOO^C  comes  out 

D  =  b  exp(-1.12  eV/kBT).  (G) 

Extrapolating  to  higher  T,  our  D  results  higher  than  that  in  [IG]  and  is  slightly  less  than  in 
[G].  All  the  data,  however,  agree  well  within  estimated  experimental  uncertainty. 

An  example  of  a  typical  response  of  the  conductivity  to  the  step  change  of  Pcd  together 
with  theoretical  curves  according  to  (5)  is  shown  in  Fig.  G.  The  logarithm  of  normal¬ 
ized  conductivity  is  plotted  there  according  to  [IG]  due  to  its  useful  linear  shape,  which  is 
conveniently  used  in  the  analysis. 


TABLE  I.  The  chemical  self- diffusion  coefficient.  Our  results  are  compared  with  data  published 
on  pure  samples  in  [6,16].  Data  with  star  were  extrapolated. 


T  (°C) 

D  (cm^/s) 

D  (cm^/s)  [6] 

D  (cm^/s)  [16] 

500 

(2.5±0.8)  X  10-^ 

*4  X  10“^ 

♦1.3  X  10"^ 

600 

(1.7±0.5)  X  10-® 

♦2.5  X  10-5 

10-5 

700 

>  4  X  10-®,  *8  X  10~® 

9  X  10-5 

5  X  10-5 

800 

♦2.5  X  10-5 

2.5  X  10-5 

1.4  X  10-5 

900 

*7.5  X  10-5 

7  X  10-5 

♦4  X  10-5 
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FIG.  6.  The  logarithm  of  the  normalized  conductivity  at  500^0  (circle)  and  600^0  (cross)  to¬ 
gether  with  theoretical  fit  (full  lines)  as  a  function  of  time  for  step-like  change  in  Pcd^  The  ctoo 
is  the  conductivity  in  an  infinite  time,  is  a  normalization  constant,  which  attributes  unity  to 
a  time,  from  which  the  curve  can  be  approximated  by  the  exponential  decay.  The  deviations  at 
short  time  are  due  to  the  surface  conduction. 

The  short-time  deviations  of  the  experiment  from  the  theory  is  due  to  the  surface  con¬ 
duction,  which  starts  to  prevail  at  low  temperatures  above  resistive  bulk.  This  effect  could 
be  disclosed  already  in  [16],  where  the  data  at  600°C  are  shown.  There  an  incorrect  inter¬ 
pretation  assuming  that  the  effect  is  a  characteristic  for  the  solution  of  diffusion  equation 
is  suggested.  As  we  see  in  Fig.  6,  the  theoretical  curve  does  not  fit  such  strong  deviation. 
Preliminary  study  shows,  that  surface  conduction  depends  on  Pcd  in  a  similar  form  as  it  is 
in  the  bulk.  The  effect  of  surface  conduction  in  medium  temperatures  has  never  been  re¬ 
ported  before.  It  complicates  the  analysis  of  experimental  data  in  the  temperature  interval 
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400“600®C,  which  is  expected  to  be  optimum  for  a  fabrication  of  low  doped  high  resistive 
CdTe.  Due  to  this  fact  the  surface  conduction  should  be  extensively  studied  in  the  future. 

Contrary  to  the  report  in  [6]  our  diffusion  data  in  increased/decreased  Pcd  produce  very 
similar  D.  This  finding  agrees  with  the  opinion  [15],  that  this  effect  is  due  to  the  capturing 
of  diffusing  species  on  doping  impurities.  In  our  pure  samples  the  interaction  with  impurities 
can  be  neglected. 

In  spite  of  experimental  divergences  in  various  runs,  we  did  not  prove  any  dependence  of 
D  on  the  stoichiometry  deviation  A.  It  implies  that  the  diffusion  coefficients  of  dominant 
point  defects  representing  excess  or  deficit  of  Cd  in  the  material  are  close  together.  This 
task  remains  open  for  future  research  in  which  the  experimental  precision  of  the  conductivity 
measurements  should  be  improved  and  the  interval  of  Pcd  far  from  P^  should  be  chosen. 

The  parallel  investigation  of  the  Hall  mobility  showed  that  the  variations  of  transport 
properties  are  mostly  operated  by  free  electrons  and  except  the  lowest  Pcd  <  O.IP^  the 
heavy  hole  transport  can  be  neglected. 

E.  Conclusion 

The  coefficient  of  chemical  self-diffusion  was  established  by  in-situ  conductivity  measure¬ 
ments.  The  data  agree  well  with  previous  results  on  pure  samples.  The  influence  of  surface 
conduction  to  the  electronic  properties  below  600^C  and  its  dependence  on  Pcd  was  found. 
This  effect  forces  an  extended  research  in  this  field  because  of  its  influence  to  the  fabrication 
of  low  doping  high  resistive  CdTe. 

IV.  CONDITIONS  FOR  FABRICATION  OF  HIGH-RESISTIVITY  CDTE 

A.  Introduction 

The  dominant  approach  for  a  production  of  high  resistivity  (HR)  detector-grade  material 
has  been  the  growth  of  Cl  doped  CdTe  from  Te-rich  solutions  mainly  by  the  THM  method. 


16 


Within  the  last  decade  substantial  effort  has  been  devoted  to  the  growth  of  high  resisitivity 
(CdZn)Te  by  the  high  pressure  Bridgman  method  (HPBM).  Here  in  some  cases  very  high 
resistivity  (10^^  ftcm)  was  achieved.  Although  HPBM  produces  large  (CdZn)Te  ingots,  their 
homogeneity  remains  an  unsolved  problem  [19].  The  manufacture  of  large  homogeneous 
ingots  has  not  reached  a  proven  stage  yet.  Continuing  effort  is  still  underway  to  improve 
the  characteristics  of  CdTe,  (CdZn)Te  materials  in  order  to  achieve  reproducible  detectors. 

One  of  the  key  problems  limiting  further  progress  in  growth  of  large,  perfect,  uniform 
and  pure  crystals  with  detector-grade  properties  is  the  insufficient  understanding  of  com¬ 
pensation  mechanisms  leading  to  the  high-resistivity  (HR)  CdTe,  (CdZn)Te.  The  demands 
on  detector-grade  material  are  comprised  in  a  low  dark  current  and  a  high  mobility-lifetime 
product  of  current  carriers.  To  prepare  such  a  material,  the  concentration  of  free  carriers 
must  be  decreased  by  a  proper  compensation  of  shallow  levels  in  the  energy  gap.  CdTe  is 
normally  synthesized  from  7N  purity  elements,  which  implies  a  level  ~  10^^  cm”^  of  foreign 
impurities.  Furthermore,  relatively  large  concentrations  of  native  point  defects  generated 
during  the  crystal  growth  can  be  frozen  in  due  to  the  fact  that  the  diffusion  effectively 
stops  at  temperature  SOO^^C  or  even  higher  depending  on  the  crystal  volume.  To  achieve 
HR  in  this  case,  the  equilibrium  of  concentrations  of  acceptors  and  donors  Na  =  Nd  with 
allowed  deviation  less  than  10^  cm”^  is  required,  which  is  difficult  to  be  reached  with  the 
state-of-the-art  technological  procedures.  Therefore  also  attempts  to  prepare  reproducibly 
HR  undoped  CdTe  controlling  the  Cd  pressure  during  the  growth  and  minimizing  this  way 
deviations  from  stoichiometric  CdTe  by  vertical  and  horizontal  Bridgman  methods  were  so 
far  not  successful!  [20] . 

Clearly,  the  pinning  of  the  Fermi  level  to  some  deep  level  near  the  middle  of  the  bandgap 
is  useful  to  achieve  the  high  resistivity  CdTe.  This  goal  is  reached  by  doping  with  impurities 
having  energy  levels  near  the  midgap  (e.g.  V,  Ti)  or  forming  complexes  with  native  defects 
at  the  midgap  levels.  Here,  however,  the  resisitivity  is  obtained  with  high  concentration 
of  defects  leading  to  the  carrier  trapping  and  recombination  and  therefore  poor  detector 
properties  of  such  samples  [21]. 
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The  purpose  of  the  present  paper  is  to  demonstrate  the  principal  difficulties  connected 
with  the  fabrication  of  HR  material  in  stoichiometric  undoped  CdTe  or  shallow  defect  doped 
CdTe.  We  shall  show,  how  the  deep  defect  doping  simplifies  this  task  and  suggest  a  techno¬ 
logical  procedure  to  prepare  high-resistivity  CdTe  with  minimized  unavoidable  deep  defect 
doping.  Vanadium  will  be  used  as  a  model  doping  element  for  demonstration. 

B.  Experiment 

Undoped  CdTe  single  crystals  with  a  diameter  up  to  100  mm  were  produced  by  the  Ver¬ 
tical  Gradient  Freeze  method  [17]  at  Cd  overpressure  1-1.3  atm.  6N  purity  starting  elements 
were  used  in  all  growth  runs.  The  cooling  speed  during  the  crystallization  process  at  tem¬ 
peratures  1120-1080'^C  was  0.5^C/hour.  All  crystals  were  very  good  quality  single  crystals 
filling  about  90%  of  the  crucible  volume  with  no  macroscopical  defects  (grain  boundaries, 
twins). 

In  situ  galvanomagnetic  measurements  of  Hall  effect  by  Van  der  Pauw  method  on  solid 
samples  with  dimensions  10x10x1.5  mm^  at  defined  Cd  pressure  Pcd  and  temperatures  up 
to  1050^C  were  performed  [3].  The  theoretical  analysis  of  such  measurements  allowed  us 
to  establish  the  high  temperature  intrinsic  conditions  in  the  P-T  diagram  and  to  correlate 
them  with  the  stoichiometry  of  the  samples. 

C.  Theory 

The  obtaining  of  high-quality  CdTe  and  (CdZn)Te  samples  for  X-  and  gamma  ray  ra¬ 
diation  detectors  is  connected  with  the  demand  to  prepare  the  HR  material  with  the  room 
temperature  (RT)  free  carrier  density  less  than  10^  cm“^,  which  is  used  as  a  limit  density 
for  the  HR  material.  The  intrinsic  concentration  at  300K  is  approximately  8  x  10^  cm~^  in 
CdTe.  The  doping  to  such  a  low  density  is  up  to  now  technologically  impossible  in  CdTe 
and  a  presence  of  any  shallow  (fully  ionized  at  RT)  uncompensated  level  results  in  the  free 
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carrier  density  well  above  the  demanded  value.  This  argument  results  in  a  necessity  to  pre¬ 
pare  the  CdTe,  where  (i)  shallow  donors  and  acceptors  compensate  mutually  {No  =  Na)  or 
(ii)  deep  defects  compensate  the  relevant  shallow  levels.  In  this  case  the  deep  donors  (accep¬ 
tors)  with  density  Ndd  {Nad)  compensate  the  uncompensated  shallow  acceptors  (donors), 
respectively.  The  Fermi  level  is  pinned  to  the  deep  level,  if  the  relations  between  densities 
Ndd  Nd  >  Na  >  Nd  or  Nad  Na>  Nd  >  Na  are  fulfilled.  It  must  be  mentioned  that 
the  identification  of  deep  defects  (donor  or  acceptor)  is  difficult  in  samples  where  shallow 
defect  density  is  unknown.  We  shall  present  the  theoretical  results  assuming  deep  donor 
doping  in  this  paper.  The  alternative  model  with  deep  acceptor  is  evident.  Simulations  of 
resistivity  in  dependence  of  the  position  of  Fermi  level  were  discussed  for  CdTe  in  Ref.  [22] 
and  for  (CdZn)Te  in  Ref.  [23]. 

The  realistic  description  of  defect  structure  in  CdTe  asks  to  include  both  extrinsic  and 
native  defects  into  the  model.  Contrary  to  the  extrinsic  defect  density  which  is  temperature 
and  pressure  independent,  the  native  defect  densities  can  be  tuned  by  a  proper  annealing 
at  defined  thermodynamic  conditions.  The  situation  can  be,  moreover,  complicated  by  an 
interaction  between  defects,  which  results  in  a  formation  of  complexes.  We  shall  neglect  the 
formation  of  complexes  first  to  simplify  the  theoretical  description.  The  effect  of  complexes 
to  the  model  will  be  discussed  subsequently. 

The  native  defects,  which  dominate  in  near  stoichiometric  samples  CdTe  are  Cd  intersti¬ 
tial  {Cdj,  divalent  donor)  and  Cd  vacancy  {Vcd,  divalent  acceptor)  [3,18].  The  equilibrium 
densities  of  neutral  native  defects  are  given  by  reaction  constants  [24] 

[cdf]  =  exp  [-^  {EiCdi)  +  UiCdir”)  +  5(C7d/)"‘']  ,  (7) 

[V^d]  =  exp  [-^  {E{Vcd)  +  U{Vcdr»)  +  ^(Vcd)”'']  .  (8) 

(9) 

Relevant  defect  formation  energy  £*,  vibrational  energy  and  entropy  in  Table  I  were 
obtained  from  the  ab  initio  calculations  and  by  fitting  the  high  temperature  transport  data 


[3], 
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rio  —  1.48  X  10^^  cm“^  is  the  number  of  unit  cells  per  volume,  is  the  Boltzmann  constant, 
and  mcd  is  the  mass  of  Cd  atom.  The  label  X  is  used  for  neutral  defects  overall. 

The  densities  of  multiply  ionized  defects  are  calculated  for  acceptors  and  donors 

[^•1  -  [^1  ^7 

where  Ea  and  Ed  are  the  acceptor  and  donor  one-electron  ionization  energies  and 
and  gxx  are  the  degeneracy  of  the  defects.  The  index  2;  means  the  ionization  degree.  The 
list  of  experimentally  determined  ionization  energies  of  native  defects  is  given  in  Table  II. 
The  experimentally  determined  energy  gap  Eg  according  to  [25]  is  used  at  all  temperatures. 
The  Fermi  energy  fip  is  obtained  solving  the  electric  neutrality  condition 

^  +  [^Cd\  +  2  [Vcd  ]  +  ~  P  [C'd/ ]  +  2  [Cd/"*"]  +  (12) 

All  extrinsic  defects  {Na,  Nd,  N^d)  are  assumed  monovalent  here.  Due  to  the  complicated 
structure  of  the  conduction  band  both  F  and  L  minima  are  taken  into  account  [10]. 

The  dynamic  of  native  defects  is  controlled  by  the  chemical  self-diffusion  described  well 
by  the  chemical  diffusion  coefficient  given  by  Eq.  6.  It  determines  the  minimum  temperature 
for  the  fast  diffusion  of  native  point  defects  in  bulk  samples.  At  low  T  the  diffusion  is  slow 
and  only  long-time  annealing  or  thin  samples  must  be  used  to  reach  the  equilibrium. 
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TABLE  11.  The  energy  and  entropy  terms  which  fit  the  high  temperature  experimental  data. 
The  underlined  parameters  were  optimized  [3].  The  others  are  from  ab  initio  calculations  [18,24]. 
Two  energies  for  the  Cdj  correspond  to  the  two  tetrahedral  interstitial  sites. 


Defect 

E  + 

S{h) 

Cdi 

0.96;  1.29 

11.1 

Vcd 

3.55 

-5.6 

TABLE  III.  The  ionization  energy  and  degeneracy  of  Cdj  and  Vcd-  Both  defects  are  divalent. 
Two  energies  for  Cdi  correspond  to  the  two  tetrahedral  interstitial  sites.  The  vacancy  energies  are 
referenced  in  the  Table,  other  data  are  according  to  [18,26]. 


Defect 

£;i(meV) 

9i 

£:2(meV) 

92 

Cdi 

0;  210 

2 

170;  360 

1 

Vcd 

50  [27] 

2 

470  [28] 
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D.  Results  and  discussion 


Before  the  start  of  the  analysis  of  doped  samples,  the  correct  description  of  undoped  CdTe 
is  necessary.  The  knowledge  of  equilibrium  concentrations  of  defects  at  high  temperature 
obtained  from  in  situ  galvanomagnetic  measurements  allowed  us  to  establish  the  position  of 
the  theoretical  stoichiometry  line  (the  so  called  p-n  line)  in  the  form  [3] 

-  8  X  10^  exp(~1.76  x  lOV^)  (atm).  (13) 

Theoretically,  the  annealing  at  should  result  in  a  fabrication  of  undoped  HR  material.  In 
reality,  however,  a  small  deviation  from  ideal  stoichiometric  annealing  conditions  results  in 
a  strong  increase  of  native  defect  density.  The  necessary  precision  of  setting  the  Cd  pressure 
at  during  the  experiment  is  practically  unrealizable  and  deep  defect  doping  is  necessary 
to  extend  the  interval  of  Pcd  in  which  HR  samples  can  be  prepared. 

The  deep  donor  level  in  the  midgap  with  the  thermal  ionization  energy  0.75  eV  below 
the  conduction  band  is  used  in  most  of  our  calculations.  Such  a  level  was  determined  in 
Vanadium  doped  CdTe  [29]  and  our  general  model  will  be  finally  demonstrated  on  samples 
doped  by  this  element. 

The  electric  neutrality  condition  (12)  is  solved  numerically  within  the  loop  of  the  P-T 
diagram  and  the  effect  of  deep  donor  doping  to  the  fabrication  of  HR  material  is  studied. 
We  show  in  Fig.  7  the  results  of  the  calculations  of  the  RT  carrier  density  in  a  deep  defect 
doped  CdTe  sample  annealed  at  600^C  and  varying  Pcd^  The  left-hand  side  of  the  figure, 
Pcd  <0.002  atm,  is  dominated  by  holes,  the  right-hand  side  shows  the  density  of  electrons. 
The  practical  impossibility  to  prepare  undoped  high- resistivity  CdTe  is  demonstrated  here. 
The  density  decreases  below  lO*^  cm“^  and  the  HR  requirements  are  fulfilled  only  in  a  very 
narrow  interval  around  the  ideal  stoichiometry  line.  The  situation  is  not  changed  also  in 
case  of  low  doping  density  =  10^^  cm“^.  The  interval  of  allowed  Pcd  broadens  and 
the  HR  conditions  are  fulfilled  more  easily  increasing  NdD‘  Concurrently,  higher  Ndd 
implicates  the  carrier  trapping  and  recombination  and  therefore  poor  detector  properties 
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of  such  samples  [21]  due  to  a  high  total  defect  density  in  the  crystal.  Thus  the  optimum 
conditions  must  be  looked  for  to  grow  the  HR  material  at  the  minimized  doping  density  with 
Pcd  being  allowed  to  spread  within  a  realistically  wide  interval.  In  view  of  the  presented 
results  the  optimum  doping  density  to  prepare  the  HR  material  by  annealing  at  600®C  gives 
Nod  in  the  region  10^^  -  10^^  cm“^.  The  plateau  in  the  hole  density  at  p  10^^  -  10^^  cm“^ 
is  due  to  the  pinning  of  the  Fermi  level  to  the  Cd-vacancy  470  meV  above  the  valence  band. 


FIG.  7.  Room  temperature  carrier  density  of  deep  defect  doped  sample  annealed  at  600^ C.  The 
black  line  corresponds  to  undoped  CdTe,  the  green,  blue  and  red  lines  correspond  to  10^^, 
and  10^^  cm“^  doping,  respectively.  The  HR  limit  density  10^^  cm“^  is  labeled  as  well. 

The  influence  of  annealing  temperature  to  the  discussed  effects  is  demonstrated  in  Fig.  8, 
where  annealing  at  500°C  was  considered.  Comparing  to  the  case  of  600°C,  the  HR  material 
can  be  really  reached  with  a  lower  deep  doping  density.  Concluding  from  that  the  annealing 
of  CdTe  at  lower  T  allows  to  reduce  the  necessary  deep  defect  doping  and  the  material  with 
better  quality  can  be  obtained.  Unfortunately,  the  annealing  at  low  T  extends  the  necessary 
time  to  reach  the  equilibrium  and  the  long-time  experiments  must  be  performed. 
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FIG.  8.  Room  temperature  carrier  density  of  deep  defect  doped  sample  annealed  at  500^ C.  The 
black  line  corresponds  to  undoped  CdTe,  the  green,  blue  and  red  lines  correspond  to  10^^, 
and  10^^  cm“^  doping,  respectively. 

The  general  impact  of  the  annealing  temperature  on  the  fabrication  of  high  resistive 
CdTe  within  the  P-T  diagram  is  demonstrated  in  Fig.  9. 

The  deep  defect  density  results  in  the  demands  to  the  thermodynamic  conditions,  which 
must  be  set  to  obtain  HR  material.  The  HR  regions  for  each  Nojo>  are  confined  by  the  full 
and  dash  boundary  lines  being  minimum  and  maximum  Pcdy  respectively.  The  labels  point 
to  the  corresponding  Ndd-  The  upper  limit  of  Pcd  for  Ndd  <  10^^  cm”^  merges  in  the 
stoichiometry  line.  In  heavily  doped  samples  =  10^^  cm~^)  HR  conditions  are  fulfilled 

overall  in  the  stability  region  below  maximum  Pcd  labeled  10^^.  We  can  deduce,  that  it 
is  possible  to  prepare  HR  samples  also  with  Ndd  approaching  to  10^"^  cm“^,  if  anneling 
temperature  T  <  450°C  is  applied. 
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FIG.  9.  Regions  in  the  P-T  diagram,  where  the  high  resistive  deep  defect  doped  CdTe  can  be 
obtained.  The  full  (dash)  lines  mark  respectively  the  minimum  (maximum)  Pcd^  which  limit  the 
HR  region.  Labels  indicate  corresponding  doping  density  in  cm“^.  For  a  plasticity,  the  HR  region 
for  doping  10^^  cm”^  is  highlighted.  The  ideal  stoichiometry  line  for  undoped  CdTe  (dash-dotted) 
is  plotted  as  well.  The  minimum  Pcd  for  doping  10^^  cm”^  results  below  the  stability  region. 

The  present  technology  does  not  allow  to  prepare  CdTe  without  shallow  extrinsic  de¬ 
fects.  Such  defects  influence  the  electric  properties,  especially  the  position  of  the  p-n  line, 
and  subsequently  the  thermodynamic  conditions  for  the  fabrication  of  HR  materials.  We 
demonstrate  the  impact  of  shallow  defects  to  the  position  of  p-n  line  in  Fig.  10,  where  its 
course  for  various  shallow  donor  or  acceptor  doping  (without  deep  defects)  is  plotted.  The 
related  RT  carrier  density  is  shown  in  Fig.  11.  Besides  the  shift  of  HR  region  we  observe 
also  its  narrowing  due  to  the  increasing  native  defect  density  far  from  P^, 
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FIG.  10.  The  effect  of  extrinsic  shallow  defect  doping  to  the  position  of  the  p-n  line  in  P-T 
diagram.  The  labels  mark  the  density  of  extrinsic  defects  for  given  p-n  line.  The  donor  doping 
shifts  p-n  line  downward  the  acceptor  doping  upward  P^.  The  dash  and  dotted  lines  show  the 
ideal  stoichiometry  line  P^  and  the  line  of  congruent  sublimation  P^,  respectively. 

The  set  of  defects  in  CdTe  should  be  completed  by  complexes,  which  are  often  mentioned 
in  the  literature.  In  our  case  a  complex  of  deep  donor  with  Vcd  should  be  taken  into  account. 
At  the  high  temperature,  at  which  the  annealing  is  performed,  the  complexes  are  probably 
unstable  and  dissociate  to  the  simple  species.  Their  involvement  is  thus  not  needed  in  the 
high  temperature  calculations.  During  the  cooling,  when  the  sample  passes  through  medium 
temperatures,  where  the  complexes  become  stable  and  the  short  range  diffusion  is  sufficiently 
fast  to  allow  the  complex  formation,  new  related  defect  levels  (possibly  shallow)  evolve  in 
the  gap.  The  complex  formation  can  have  two  inevitable  consequences  for  the  HR  samples. 

(i)  The  density  of  single  deep  donors  decreases  as  they  are  lost  by  the  complex  formation. 

(ii)  New  complex  levels  act  as  new  shallow  acceptors  or  donors  and  must  be  included  into 
the  electric  neutrality  calculation  (12).  Consequently,  the  Fermi  level  could  drop  off  the 
midgap  and  the  resistivity  can  significantly  decrease. 
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FIG.  11.  Room  temperature  carrier  density  of  sample  annealed  at  500^C  for  various  shallow 
doping  density  at  Ndd  =  5  x  10^^  cm”^.  The  dash  line  plots  the  zero  shallow  doping,  the  labeled 
lines  a,  b,  c  correspond  to  10^^,  10^®,  10^^  cm"^  shallow  donor  doping,  d,  e,  f  to  10^^,  10^®,  10^"^ 
cm~^  shallow  acceptor  doping,  respectively. 

For  a  simple  estimation  within  our  model,  we  assume,  that  shallow  acceptor  complex  X 
comes  up  via  reaction  of  deep  donor  with  Vcd  ^nd  moreover  all  Vcd  enter  in  the  complex  or 
recombine  with  Cdi,  i.e.  [X]  =  [Vcd]  -  [Cd^],  and  foreign  shallow  defects  being  neglected. 
Then  the  criteria  for  HR  sample  without  complex  Nj^d  >  2([Vcd]  —  [Cdj])  >  0  transform 
into  Ndd  —  [X]  >  z[X],  where  2:  is  the  valence  of  the  complex.  If  2:  =  1,  the  situation  is 
not  influenced  by  the  complex  formation  and  the  sample  remains  high  resistive.  If  z  >  1, 
the  region  for  the  fabrication  of  HR  material  in  P-T  diagram  (Fig.  9)  narrows  and  the  deep 
defect  density  must  be  increased  or  annealing  temperature  decreased  to  fit  the  HR  region 
well. 

Vanadium  is  often  used  as  a  typical  deep  defect  in  CdTe.  The  position  of  its  level  formally 
identified  to  the  single  donor  in  Cd  site  was  found  in  the  midgap  with  the  thermal 

ionization  energy  0.75  eV  below  the  conduction  band  [29].  If  such  an  interpretation  is  correct, 
two  shallow  donor  levels  of  V  must  be  compensated  by  native  or  foreign  acceptors  in  Cd  sites 
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to  prepare  HR  V-doped  material.  In  view  of  the  previous  discussion  the  compensation  by 
foreign  shallow  acceptors  is  preferred.  They  compensate  the  material  near  the  stoichiometry 
line  where  the  tuning  of  the  position  of  the  Fermi  level  by  the  variation  of  Pcd  is  the  most 
effective.  For  the  practical  demonstration  we  show  in  Fig.  12  the  resistivity  of  V-doped 
sample  annealed  at  500®C  with  ideally  compensated  shallow  levels. 


FIG.  12.  The  room  temperature  resistivity  of  V-doped  sample  with  ideally  compensated  shallow 
levels  annealed  at  500^ C  as  a  function  of  annealing  Cd  pressure.  The  red  and  blue  lines  correspond 
to  the  V-doping  densities  10^®  and  10^^  cm~^,  respectively. 


E.  Conclusion 

Methods  for  fabrication  of  high  resistivity  detector-grade  CdTe  are  discussed.  Thermo¬ 
dynamic  conditions  at  which  an  undoped  material  is  stoichiometric  are  found  based  on  high 
temperature  in  situ  galvanomagnetic  measurements,  when  two  major  native  defects  -  the 
Cd  vacancy  and  the  Cd  interstitial  are  applied  in  model  calculations.  It  is  demonstrated, 
that  only  a  small  deviation  from  stoichiometric  conditions  results  in  a  strong  generation  of 
native  defects  and  therefore  in  a  crystal  resitivity  decrease.  This  makes  the  fabrication  of  a 
high  resistivity  undoped  CdTe  a  problem,  which  is  difficult  to  solve  with  the  state-of-the-art 
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technological  possibilities.  Modelling  of  room  temperature  crystal  resistivity  in  dependence 
of  the  concentration  of  a  dopant  with  the  energy  level  near  the  middle  of  the  gap  and  thermo¬ 
dynamic  conditions,  from  which  the  crystal  is  quenched  to  room  temperature  is  performed. 
It  is  shown,  that  when  increasing  the  deep  defect  density  the  range  of  optimal  Pcd  at  each 
annealing  temperature  broadens  and  the  high  resistivity  can  be  reached  more  easily.  The 
influence  of  shallow  extrinsic  defects  and  complexes  of  native  defects  with  the  deep  dopant 
results  in  shifting  and  narrowing  of  the  Pcd{T)  range,  where  the  high  resistivity  material 
can  be  obtained.  The  presented  model  can  be  generalized  to  any  structure  of  defect  levels 
in  CdTe. 


V.  DEFECT  INDUCED  OPTICAL  TRANSITIONS  IN  CDTE 

While  at  high  temperatures  in  CdTe  native  point  defects  play  a  dominant  role  in  not 
intentionally  doped  crystals,  room  and  low  temperature  optical  and  galvanomagnetic  proper¬ 
ties  are  often  strongly  influenced  by  residual  impurities.  The  photoluminescence  is  a  favorite 
method  used  for  a  study  of  defect  structure.  Systematic  effort  was  devoted  to  And  out  the 
binding  energy  of  shallow  acceptors  by  low  temperature  photoluminescence  (PL),  PL  ex¬ 
citation  spectroscopy  and  magnetooptical  analysis  [30-32].  Doping  experiments  performed 
on  high  purity  materials  have  shown,  that  all  studied  acceptors  are  related  to  substitutional 
impurities  [30].  Contrary  to  a  large  amount  of  papers  devoted  to  photoluminescence  on 
CdTe  and  related  compounds,  optical  absorption  studies  are  rare.  Infrared  absorption  (400- 
2000  cm”^)  in  Cu  doped  CdTe  was  analyzed  in  [33],  while  combined  photoluminescence  and 
absorption  measurements  of  the  1.1  eV  deep  band  related  to  Te  vacancies  were  performed 
in  [34].  Irregularities  in  the  structure  of  absorption  edge  in  CdTe  were  explained  in  [35] 
by  a  phonon  assisted  direct  excitonic  absorption.  The  shape  of  the  absorption  spectra  as 
well  as  the  interpretation,  which  agrees  qualitatively  with  our  paper,  was  reported  for  ZnTe 
in  [36].  In  this  paper  we  evaluate  the  absorption  below  the  absorption  edge  of  CdTe  and 
Cdi_a;Zna:Te.  We  use  a  model  of  totally  charged  classical  acceptor  impurity  band,  where  the 
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Coulomb  interaction  is  screened  by  mobile  charges  in  partially  filled  level. 

A.  Experimental 

CdTe  and  CdZnTe  (with  the  Zn  content  ^^4%)  single  crystals  were  fabricated  in  our  lab¬ 
oratory  by  the  Vertical  Gradient  Freeze  method  [17].  6N  purity  starting  elements  were  used 
in  all  growth  runs.  The  Cd  pressure  at  the  maximum  temperature  during  the  growth  was 
reached  by  adding  of  the  corresponding  amount  of  Cd  in  the  ampoule.  All  crystals  were  very 
good  quality  single  crystals  filling  about  90%  of  the  crucible  volume  with  no  macroscopical 
defects  (grain  boundaries,  twins),  p-type  with  the  room  temperature  concentration  ^  10^^ 
cm"^.  No  intentional  doping  was  applied.  For  our  PL  and  optical  absorption  studies  samples 
with  thickness  1.5-3  mm  were  prepared  by  cutting,  polishing  and  etching  in  Br-methanol 
solution. 

The  transmittance  spectra  (resolution  1  cm”^)  were  measured  by  Fourier  transform  spec¬ 
trometer  IFS66/s  (Bruker)  and  He  flow  cryostats  VSK  4-300  (Leybold). 

B.  Theory 

The  structure  of  the  absorption  edge  in  semiconductors  is  generally  determined  by  four 
effects  which  are  responsible  for  the  finite  absorption  constant  a  below  the  principal  band- 
edge  energy  Eg.  1.  The  exciton  absorption  is  observed  below  Eg  forming  free  exciton  with 
binding  energy  in  CdTe  10.6  meV,  and  the  donor  and  acceptor  bound  exciton  at  13  meV 
and  16.8  meV,  respectively.  2.  The  band  edge  is  broadened  due  to  the  lattice  disorder, 
which  results  in  the  exponential  tail  in  the  density  of  states  and  is  manifested  also  in  a. 
3.  The  lattice  defects  initiate  localized  states  in  the  gap.  The  absorption  connected  with 
the  electron  transfer  into/from  these  levels  allows  to  characterize  them  and  to  establish 
their  ionization  energies.  4.  All  transitions  can  be  accompanied  by  the  phonon  assisted 
processes,  at  low  temperature  characterized  by  the  emission  of  optical  phonon  at  the  energy 
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21.2  meV  in  CdTe.  The  effect  results  in  the  repeated  characteristic  structure  accompanying 
the  transitions  in  item  3. 

The  purpose  of  this  paper  is  to  explain  the  absorption  20  -  70  meV  below  Eg  at  low 
temperature  (4~20K),  where  the  remarkable  details  were  found  (see  Fig.  13).  In  respect 
to  the  items  1.-4.  above  we  can  eliminate  the  excitonic  effects,  which  are  not  observed  at 
such  energies.  The  exponential  broadening  of  the  absorption  edge  forms  the  background 
absorption  in  this  region.  The  characteristic  wavy  shape  is  not,  however,  explained  by  this 
mechanism.  Based  on  given  remarks,  we  suggest  that  the  mechanism  in  item  3  can  be 
expected  to  describe  the  observed  effect. 

The  single  impurity<-^band  absorption  constant  is  calculated  by  a  standard  way  according 


to  [37,38]  in  the  form 


ai{hcj)  — 


3ceQnh  CO  \2m*a^ 


yhco  —  Eg  -f-  Ej 
— h  hco  —  Eg  -{-  Ej 


where  P  =  8.5  x  10"®  eVcm  is  the  momentum  matrix  element  [39],  n  is  refractive  index, 
and  e,  c,  ^o,  h  are  elemental  constants.  We  shall  show  in  the  discussion  that  the  absorption 
is  connected  with  the  acceptor-conduction  band  (CB)  transitions,  where  the  acceptors  are 
charged  due  to  compensating  donors.  In  this  case  the  acceptor  degeneracy  g  —  A,  Ej  (=  Ea) 
and  as  are  the  acceptor  binding  energy  and  Bohr  radius  in  the  effective  mass  approximation 
and  m*  —  0.096mo  is  the  CB  effective  mass. 

At  the  low  temperature,  when  the  free  particle  density  can  be  neglected,  the  potential 
at  the  localized  states  is  influenced  by  the  Coulomb  interaction  with  charge  on  other  point 
defects.  Consequently,  Ea  fluctuates  around  its  ideal  single- acceptor  value  and  the  impurity 
band  (IB)  with  the  density  of  states  G{E)  is  formed.  The  energy  deviation  of  the  respective 
acceptor  is  obtained  in  the  form 

dEA{v)  =  (15) 

ATTSoSr  |r  —  Tsl 

where  Zs  is  the  charge  of  the  respective  site.  The  total  absorption  constant  a  is  obtained 
integrating 
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(16) 
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/  aiihw  +  X  -EA)GiX)dX, 

J  Eg— hw 

where  G{X)  influences  the  start  of  the  absorption  which  is  smoothed  comparing  to  the  sharp 
increase  in  case  of  isolated  acceptor  level.  The  fitting  of  a  to  experimental  data  allows  to 
estimate  the  shallow  defect  parameters.  The  magnitude  of  a  determines  the  defect  density 
and  the  start  of  the  absorption  gives  the  position  of  respective  ionization  level.  The  shape 
of  the  absorption  edge  is  used  for  the  examination  of  the  relevant  IB. 

The  impurity  levels  in  CdTe  are  of  very  complicated  structure  forming  both  shallow  and 
deep  donor  and  acceptor  states.  The  exact  distribution  of  the  levels  is  not  known  in  our 
samples.  In  respect  to  this  fact  we  use  a  minimal  model  here,  which  includes  a  set  of  shallow 
acceptors  to  describe  the  characteristics  in  the  absorption,  and  one  compensating  donor 
level  (does  not  matter  if  deep  or  shallow),  placed  sufficiently  far  above  the  valence  band. 
While  acceptors  are  assumed  ionized  (-),  donors  are  partly  neutral  to  fulfill  the  neutrality 
condition  and  the  Fermi  energy  is  placed  near  this  level.  The  charge  on  donors  (+)  occupies 
preferentially  the  sites  near  acceptors  and  screens  the  Coulomb  interaction.  Its  equilibrium 
distribution  in  finite  temperature  is  obtained  by  Monte  Carlo  simulations  [40-42], 

The  final  states  in  CB  are  also  disturbed  by  potential  fluctuations.  Their  description 
is,  however,  much  more  complicated  and  runs  to  the  necessity  to  solve  the  problem  of 
localization  in  disordered  systems.  The  solution  of  this  task  is  out  of  the  scope  of  this 
paper.  We  estimate  such  effect  here  using  the  well-known  relation  between  binding,  kinetic 
and  potential  energy  Eb,  <  T  >  and  <  V  >,  respectively,  which  reads  in  the  hydrogen¬ 
like  states  Eb  =<  T  >=  —  <  V  >  /2.  Relying  to  that  we  assume  that  the  potential 
fluctuations  of  the  CB  edge  manifest  by  one  half  in  the  transition  energy  and  decrease  twice 
the  effective  width  of  G{E)  in  comparison  to  Eq.  (15).  On  the  contrary,  the  opposite  action 
to  G  can  be  observed  due  to  the  unknown  charged  defects,  both  donors  and  acceptors, 
which  increase  the  potential  fluctuations  and  result  in  the  increase  of  the  width  of  G.  The 
shape  of  the  absorption  edge  can  give  an  information  about  G,  which  can  be  thus  used  for 
the  characterization  of  the  samples  with  respect  to  their  purity  and  make  visible  by  the 
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indirect  way  the  charged  defects,  which  are  not  observed  in  other  experimental  methods. 
In  our  calculations  we  neglect  first  these  effects  assuming  flat  CB  and  include  only  the 
acceptor  levels,  which  are  observed  in  the  absorption.  The  above  mentioned  corrections  are 
subsequently  discussed  to  become  the  more  realistic  model. 

C.  Results  and  discussion 

The  first  task,  which  is  to  be  solved  is  the  decision  about  the  type  of  the  point  defects 
which  are  responsible  for  the  sub-gap  transitions.  As  it  is  seen  in  Fig.  13,  the  absorption 
can  be  good  approximated  by  the  square  root  of  energy  in  agreement  with  the  numerator  in 
Eq.  (14),  The  denominator  must  have  only  weak  influence.  Such  demand  is  easy  fulfilled  in 
case  of  acceptors  for  which  a%)  ~  0.3  eV  for  hydrogen-like  acceptor  is  significantly 

greater  than  the  energy  interval  of  our  interest.  The  opposite  situation  is  met  in  case  of 
transitions  from  valence  band  to  donors,  for  which  the  factor  above  is  significantly  suppressed 
and  a  would  have  to  pass  maximum  within  the  monitored  interval.  Another  argument  can 
be  based  on  experimental  findings  [43],  which  report  acceptor  levels  close  to  our  observations. 
No  donor  level  was  reported  in  this  region. 
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FIG.  13.  Absorption  and  luminescence  spectra  of  (CdZn)Te  at  4.8K.  The  theoretical  a  is  not 
resolved  from  experiment  in  this  scale. 

We  present  results  of  the  fit  of  a  on  two  our  samples  measured  at  4.2  K.  A  very  good 
agreement  of  theoretical  fit  with  experimental  absorption  together  with  luminescence  of 
sample  Cdo.96Zno.04Te  is  presented  in  Fig.  13,  for  which  Eg  =  1.625  eV  was  established. 
We  see  the  characteristic  advancement  of  a  near  1.57  eV,  which  is  greatly  fitted  within 
presented  model.  The  position  of  level  Eai  —  58.6  meV  agrees  well  with  Eas  =  58  ±  2 
meV  and  Esb  =  61  ±  2  meV  [43].  The  interpretation  is  also  supported  by  the  position  of  the 
donor-acceptor  peak  DAP  in  the  luminescence,  which  is  shifted  to  lower  energy  by  the  donor 
binding  energy.  The  acceptor  density  at  this  level  Nai  =  2.8  x  10^^  cm“^  was  calculated 
according  to  (14),  (16),  where  the  Bohr  radius  aui  =  10. 5A  was  used  [43]. 

The  similar  motif  is  repeated  once  more  at  higher  energy  1.59  eV.  We  tried  to  describe 
this  part  by  a  new  acceptor  level  with  Ea2  =  38.2  meV,  Such  level  is  not  reported  in  the 
literature  and  its  Bohr  radius  aB2  is  thus  not  known.  We  connected  aB2  with  ajgi  by  the 
relation  aB2  =  clb\\J Eai/ Ea2,  which  yields  aB2  =  13.0A  and  resulting  acceptor  density 
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Na2  =  3.2  X  10^^  cm-^ 

In  the  luminescence  spectrum  the  standard  free  exciton  (X),  donor-  (D,X),  acceptor- 
bound  excitons  (A,X)  and  their  phonon  replicas  are  observed.  A  broad  luminescence  band 
at  1.55-1. 58eV  includes  the  DAP  peak,  the  second  phonon  replica  of  the  (A,X)  and  the  peak 
corresponding  to  a  transition  from  acceptor  to  the  conduction  band  (e,A).  Due  to  the  fact, 
that  these  contributions  to  the  luminescence  band  are  not  well  resolved,  determination  of 
the  energy  level  of  the  acceptor  from  the  luminescence  spectrum  is  difficult. 

The  detailed  analysis  of  the  results  of  the  Monte-Carlo  simulations  on  this  sample  is 
shown  in  Fig.  14,  where  the  theoretical  a,  background  absorption  and  the  density  of  states 
G  of  the  impurity  bands  is  presented.  We  fit  the  background  a  by  a  linear  and  exponential 
functions,  which  fit  properly  peripheral  parts  of  the  investigated  interval.  The  course  of  G 
of  both  levels  is  the  same  due  to  the  same  character  of  Coulomb  interaction  in  our  classical 
approach  with  the  line  width  of  FWHM=5.6  meV.  The  best  fit  was  obtained  with  the 
compensating  donor  density  Nd  =  9  x  10^"^  cm“^,  whose  variations  influence  slightly  the 
shape  of  G. 

The  excellent  fit  was  obtained  within  our  very  simple  model  of  two  shallow  acceptors 
and  one  compensating  donor.  If  the  additional  charges  are  included  to  the  model,  the  IB 
broadens  and  the  start  of  the  absorption  is  smoothed.  The  influence  of  charged  levels  to  a  is 
demonstrated  in  Fig,  15,  where  new  fully  charged  donor  and  acceptor  levels  with  the  same 
density  were  supplemented  in  the  model  calculations  and  the  compensation  of  the  sample 
was  increased.  The  finding  implies  that  the  sample  is  pure  without  substantial  concentration 
of  other  invisible  charged  defects.  The  situation  could  be  complicated  by  the  warping  of  the 
conduction  band  discussed  in  the  Theory.  If  G  is  narrowed  twice,  the  best  fit  is  obtained 
with  the  invisible  charge  density  ~  4  x  10^^  cm~^.  The  careful  study  of  CB  in  disordered 
systems  is  thus  desirable. 
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Energy(eV) 

FIG.  14.  Theoretical  absorption  of  the  (CdZn)Te  sample  (full  line  -  1),  background  absorption 
(dash  '  2)  and  the  absorption  connected  with  acceptor-conduction  band  transition  (dash-dotted  - 
3).  The  line  1  is  the  sum  of  lines  2  and  3.  The  peaks  show  the  impurity  bands  density  of  states 
calculated  by  Monte-Carlo  simulations  on  a  cubic  cell  containing  4000  acceptors. 
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FIG.  15.  The  influence  of  the  background  doping  to  the  shape  of  the  absorption.  The  legend 
indicates  the  concentration  of  the  additional  donor  Nd  and  acceptor  Na  charged  levels  {Nd  =  Na)- 

The  interpretation  of  the  second  level  Ea2  is,  however,  not  clear  and  alternative  descrip¬ 
tions  could  be  also  valid  in  this  case.  The  corresponding  donor-acceptor  peak  is  not  observed 
in  the  luminescence  and  the  ionization  energy  difference  Eai  —  Ea2  ~  20.4  meV  is  close  to 
the  optical  phonon  energy  21.2  meV.  This  finding  points  to  a  possibility  to  interpret  this 
effect  as  a  phonon  assisted  transition  from  the  acceptor  Eai^  Other  plausible  explanation 
could  be  accepted  connecting  the  effect  with  the  complicated  structure  of  the  band  edge, 
which  is  responsible  for  the  strong  increase  of  a  above  1.60  eV. 

The  results  of  the  measurements  and  the  fit  of  a  on  sample  CdTe  are  plotted  in  Fig.  16 
{Eg  —  1.605  eV).  The  typical  character  of  the  absorption  edge  is  also  evident  at  1.55  eV 
and  the  fit  produces  Eai  —  59-0  meV,  again  in  good  relation  with  DAP.  The  level  density 
results  Na\  =  1-4  x  10^^  cm”^  in  this  case  and  the  density  of  states  FWHM=2.8  meV.  Other 
acceptor  level  corresponding  to  the  absorption  edge  near  1.57  eV  is  not  well  pronounced 
and  cannot  be  fit  within  our  model.  As  in  the  previous  sample  its  position  is  again  about 
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20  meV  above  £*^11,  which  supports  the  interpretation  of  this  part  by  the  mechanism  of  the 
phonon  assisted  absorption  from  the  level  £ai- 


FIG.  16.  Absorption  and  luminescence  spectra  of  CdTe  at  4.8K.  The  theoretical  a  deviates  from 
experiment  at  1.565  eV  due  to  the  complicated  course  of  the  absorption  edge,  which  is  not  included 
in  the  model. 

Comparing  shape  of  the  fundamental  absorption  edge  in  both  samples,  we  observe  the 
start  of  strong  increase  of  o;/  at  similar  energies  relatively  to  Eg,  In  CdTe  the  fundamental 
absorption  starts  at  about  20  meV  lower  energy  relatively  to  Eg  and  corrupts  partly  the  shape 
of  a  in  the  region  of  defect  induced  absorption.  A  seeming  contradiction  between  FWHM  of 
G,  which  implies  a  higher  purity  of  the  CdTe  sample  and  the  opposite  result  obtained  from 
the  shape  of  the  absorption  edge  could  be  explained  by  a  higher  concentration  of  neutral 
defects  in  CdTe  sample.  Neutral  defects  do  not  manifest  in  the  Coulomb  broadening  of  G, 
they  can,  however,  influence  the  band  edge.  The  width  of  the  impurity  band  enables  to 
estimate  the  total  concentration  of  charged  impurities  in  the  sample  including  deep  defects, 
which  are  difficult  to  be  observed  by  other  analytical  methods.  The  concentration  of  charged 
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defects  in  both  samples  is  much  lower  than  it  could  be  deduced  from  the  6N  purity  starting 
elements.  This  result  concludes,  that  foreign  impurities  are  mostly  neutral  or  form  neutral 
complexes  in  our  samples. 

VI.  REFRACTION  INDEX  OF  CDTE:  SPECTRAL  AND  TEMPERATURE 

DEPENDENCE 

CdTe  has  considerable  application  in  radiation  detectors  and  other  optoelectronics  de¬ 
vices  and  is  used  as  substrate  for  epitaxial  growth,  e.g.  of  HgCdTe  layers  for  infrared  de¬ 
tectors  fabrication.  Therefore,  properties  of  CdTe  have  been  investigated  intensively.  Many 
papers  (e.g.  [35,44-48])  reported  on  a  spectral  dependence  of  refractive  index  at  room  tem¬ 
perature  (RT),  a  lot  of  works  dealt  with  optical  parameters  above  the  gap  energy  at  RT 
as  well.  Temperature  dependence  of  the  energy  gap  Eg  has  been  investigated  by  various 
methods  (  [49]  and  references  therein).  Temperature  dependence  of  the  spectral  position  of 
a  higher  critical  point  Ex  was  measured  by  Ksendzov  et  al  [50].  Optical  parameters  above 
the  energy  gap  were  investigated  at  elevated  temperatures  (from  RT  to  727  K)  by  Kim  and 
Sivananthan  [51].  To  our  knowledge,  no  systematic  study  of  temperature  dependence  of 
refractive  index  has  been  published. 

The  data  reported  in  this  paper  were  obtained  by  measurement  of  interference  oscillations 
in  transmittance  for  20  temperatures  between  10  K  and  295  K.  We  tried  to  describe  both 
the  spectral  and  temperature  dependence  of  the  refractive  index  in  the  transparency  range 
by  simple  empirical  relations. 

The  sample  used  in  our  study  was  nominally  undoped,  high  resistivity  bulk  CdTe.  It 
was  ground  and  polished  to  the  thickness  54.6  //m,  as  determined  by  the  surfometer  SF 
200  (Planer  Products  Ltd.).  The  transmittance  spectra  (resolution  1  cm”\  0.12meV)  were 
measured  using  Fourier  transform  spectrometer  IFS66/s  (Bruker)  and  helium  flow  cryostats 
VSK  4-300  (Leybold).  The  divergence  of  the  beam  was  decreased  by  a  diaphragm  down  to 
4®.  At  temperature  lOK,  430  interference  maxima  were  recognized  in  the  investigated 
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FIG.  17.  Spectra  of  refractive  index  for  temperatures  lOK,  150K  and  295K. 

spectral  range,  while  377  maxima  were  observed  at  295K  due  to  a  narrower  spectral  range 
of  the  sample  transparency.  The  transmittance  spectra  near  the  absorption  edge  had  been 
treated  by  a  baseline  procedure  before  the  positions  of  the  interference  maxima  were  estab¬ 
lished.  Each  maximum  was  fitted  by  a  cosine  -  like  function  to  determine  its  position  more 
accurately.  Indices  of  refraction  were  calculated  by  well  known  formula  for  a  position  of 
interference  maxima.  The  obtained  spectra  for  three  temperatures  are  shown  in  Fig.  17. 

We  compared  our  experimental  results  at  RT  with  some  of  published  relations.  The 
differences  reach  several  hundredths  in  the  spectral  region  near  the  absorption  edge,  see 
Fig.  18a, b.  Marple  [35],  Wei  et  al.  [47]  and  Peiris  et  al.  [48]  used  Sellmeier’s  formula  for 
dispersion  caused  by  one  oscillator.  The  respective  relations  were  obtained  by  fitting  of 
the  experimental  values  in  various  spectral  ranges:  Marple  [35]  0.5  -  1.4  eV  by  method  of 
minimum  deviation  at  dispersion  by  the  prism;  Wei  et  al  [47]  0.75  -  1.5  eV  from  ellipsometric 
data;  Peiris  et  al.  [48]  1.2-1.45  eV  from  oscillation  of  reflectivity  on  a  MBE  epitaxial  layer 
and  by  a  prism  coupler  technique  at  fixed  wavelength  852.3  nm. 
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(a) 


FIG.  18.  (a),(b)  Differences  between  our  experimental  data  and  values  calculated  from  published 
empirical  relations,  (c)  Crosses  represent  residuals  in  fitting  procedure  of  our  data  at  295K  to 
Sellmeier’s  formula  (1-oscillator  model).  Solid  line  shows  differences  between  the  experimental 
data  and  calculation  using  Eqs.  (17-22). 


Finkman  et  al  [44]  supplied  Sellmeier’s  formula  by  a  term  describing  a  decrease  of  the  index 
at  low  photon  energies  due  to  phonons.  Pikhtin  and  Yaskov  [46]  and  Kucera  [45]  used  in 
their  fitting  equations  terms  derived  by  the  Kramers  -  Kronig  transformation  in  a  model  of 
rectangular  spectrum  of  imaginary  part  of  dielectric  function  to  take  into  account  the  band 
character  of  the  energy  spectrum. 

The  results  of  fitting  of  our  experimental  results  by  Sellmeier’s  formula  can  be  seen  in 
Fig.  18c  for  RT.  A  better  agreement  with  experiment  can  be  obtained  by  increasing  the 
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number  of  the  oscillators.  It  has  been  shown  that  model  taking  4  oscillators  into  account 
(hereafter  model  S4)  describes  the  spectra  quite  well.  There  is  a  lot  of  parameters  in  the 
model,  some  of  them  are  interdependent  in  fitting  procedures.  We  tried  various  assumptions 
concerning  both  the  temperature  dependence  of  oscillators  energies  and  their  strengths. 
Below  we  present  relations  obtained  from  a  model  (S4),  where  the  oscillators  strengths  were 
assumed  to  be  constant,  excepting  the  strength  of  the  ’’exciton”  oscillator.  Furthermore,  the 
spectrally  independent  term,  the  energy  and  strength  of  the  phonon  oscillator  could  be  fixed 
as  well.  After  this  simplification  the  form  with  20  parameters  was  fit  by  the  least  square 
method  to  the  experimental  data 

where  a  =  3.153,  the  term  z  =  1  denotes  the  contribution  of  optical  transitions  near  band 


=  1.5976-5^1^:1^,  lev] 


9i{T)=  0.0369- 


9.26  X  10-®r^ 
T  +  22.4 


,  0.0946  \  . 

^  l  +  exp(20-^)j’f®  ^ 


terms  with  i=2  and  i=3  approximate  an  effect  of  all  higher  transitions 


=  1.7102 
£3(r)  =  3.2668 


.92  =  0.3454  eV^ 
93  =  39.966  eV^ 


and  term  i  =  A  represents  lattice  vibrations 


Ei  =  0.018  eV 


=  9.02  X  10"^  eV^ 


The  precision  of  the  fitted  parameters  is  given  with  respect  to  their  influence  to  the  least 
square  sum.  It  is  surprising  that  this  simplified  model  describes  the  experimental  data  so  well 
with  a  maximum  deviation  near  the  absorption  edge  (about  20  meV  below  Eg)  An  —  ±0.002. 
The  error  of  n  due  to  an  uncertainty  of  the  thickness  of  the  sample  is  estimated  to  ±0.5%. 
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The  rise  of  the  index  with  increasing  temperature  is  caused  by  a  decrease  of  the  transition 
energies  in  a  qualitative  agreement  with  published  dependence  of  gap  energy  [49]  and  energy 
of  higher  critical  points  [50,51],  The  shifts  of  our  parameters  E2  and  Es  are  weaker  than  the 
reported  shifts  of  the  critical  points. 

Temperature  dependence  of  Ei  is  slightly  different  from  the  exciton  energy  (determined 
from  the  luminescence  spectra)  as  shown  in  Fig.  19a.  We  let  the  parameter  gi  free  during  the 
fitting  procedure  to  obtain  a  good  agreement  with  experimental  data.  The  non-monotonic 
temperature  dependence  with  a  well  pronounced  step  near  T=30K  shown  in  Fig.  19b  was 
approximated  by  a  Fermi-Dirac-like  term  in  Eq.  (19).  Irregularities  near  30K  are  observed 
in  the  absorption  as  well.  Spectral  position  of  the  absorption  edge  defined  by  the  absorption 
coefficient  500  cm  ^  is  shown  in  Fig.  19a.  These  effects  are  caused  by  a  temperature 
dependence  of  state  population  near  the  band  edges  due  to  irradiation  of  the  sample  in 
the  spectrometer.  We  have  measured  the  position  of  the  absorption  edge  at  extremely  low 
intensity,  which  resulted  in  a  shift  of  the  edge  down  by  2  meV  at  lOK,  see  dashed  line 
in  Fig.  19a.  We  suppose  that  this  correction  could  modify  the  fitted  gi  by  removing  the 
Fermi-Dirac-like  term  in  (19).  The  exact  verification  of  such  model  is,  however,  very  difficult 
in  our  experimental  setup  because  of  noise  in  the  interference  oscillations  at  low  radiation 
intensity. 
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FIG.  19.  (a)  Parameter  Ei{T)  compared  with:  exciton  energy  as  determined  from  luminescence 
of  our  sample;  with  the  relation  Eg{T)  given  by  [49];  with  the  absorption  edge  defined  by  the  value 
500  cm“^  determined  from  measurement  at  usual  intensity  of  radiation  in  the  spectrometer  (big 
circles)  and  at  intensity  1000 x  lower  (small  circles),  (b)  Temperature  dependence  of  the  parameter 
gi  as  given  by  Eq.  (19)  (solid  line)  and  gi  expressed  by  the  first  bracket  (without  Fermi-Dirac 
term)  in  Eq.(19)  (dashed). 

The  temperature  dependences  of  the  refractive  index  are  shown  in  Fig.  20  together  with 
the  curves  calculated  from  equations  (17)- (22)  for  several  wavelengths.  The  experimental 
values  of  refractive  index  at  these  wavelengths  were  determined  by  an  interpolation  between 
the  nearest  interference  maxima. 
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FIG.  20.  Temperature  dependence  of  refractive  index  for  indicated  wavelengths.  Solid  lines  are 
calculated  by  Eqs.  (17)- (22). 

We  can  summarize  that  simple  empirical  relation  can  be  used,  approximating  spectral 
and  temperature  dependence  of  the  refractive  index  of  CdTe  with  a  good  accuracy  with 
maximum  deviation  An  =  ±0.002  near  the  absorption  edge.  This  form  can  be  used  in 
practical  applications  and  as  a  check  in  formula  for  a  formulation  of  more  realistic  model 
complying  with  the  band  structure  of  CdTe. 
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